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The phase diagram of the two-leg t-Jz ladder is explored, using the density matrix renormaliza- 
tion group method. Results are obtained for energy gaps, electron density profiles and correlation 
functions for the half-filled and quarter-filled cases. The effective Lagrangian velocity parameter Vp 
is shown to vanish at half-filling. The behaviour of the one-hole gap in the Nagaoka limit is investi- 
gated, and found to disagree with theoretical predictions. A tentative phase diagram is presented, 
which is quite similar to the full t-J ladder, but scaled up by a factor of about two in coupling. 
Near half-filling a Luther-Emery phase is found, which may be expected to show superconducting 
correlations, while near quarter-filling the system appears to be in a Tomonaga-Luttinger phase. 
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I. INTRODUCTION 

The discovery of high-temperature superconductiv- 
ity in the cuprate materials has sparked huge inter- 
est in models of strongly correlated electrons in low- 
dimensional systems, such as the Hubbard, t-J and t- 
Jz models. These models are exactly solvable in one 
dimension, at least in some special cases; but the two- 
dimensional models pose a formidable numerical chal- 
lenge. The 'minus sign' problem is a major stumbling 
block for Monte Carlo calculations in these fermionic sys- 
tems; and the convergence of Density Matrix Renormal- 
ization Group (DMRG) calculations is slow in two di- 
mensions. Exact finite-lattice calculations are limited to 
small lattice sizes; while series expansions have typically 
only been useful for special cases such as the half-filled 
limit. 

In these circumstances, a considerable effort has been 
invested in the study of 'ladder' systems consisting of 
two or more coupled chains of sites. Ladders provide 
a 'half-way house', in some sense, between one and two 
dimensions; and they also display some very interesting 
effects in their own right They display quite differ- 

ent behaviour depending on whether the number of legs 
is even or odd, as in the Haldane effect for the Heisenberg 
ladders Furthermore, experimental compounds have 
been found which form ladders such as SrCu203 0, 
which may allow the theoretical predictions to be tested 
exp er iment ally. 

The t-J model is an 'effective Hamiltonian' for the par- 
ent Hubbard model, valid when the Coulomb repulsion 
is large, but nowadays it is considered as an interesting 
model in its own right The t-Jz model is a vari- 

ant in which the rotational symmetry is broken, and the 
spin interactions are Ising-like. The two-leg t-J ladder 
has been extensively studied, using exact diagonalization 
millilll, 12, 13, 14,15], quantum Monte Carlo [H, 
the DMRG technique 111 (llMM III 111 El , a com- 
bination of different methods '24l|, or using mean-field 



or approximate analytic methods 0, |2^ |26L . Near 
half-filling, the model has been explored using dimer se- 
ries expansions |2^ IH |33| . 

Our object is to study the corresponding two-leg t- 
Jz ladder, and compare the results for the two models. 
This model has not been studied before, as far as we are 
aware, but the t-Jz chain has been discussed by Batista 
and Ortiz |^, and the t-Jz model on the square lattice 
has been treated by several groups IH lU 0, 113 . Our 
primary tool is the DMRG approach, supplemented with 
a few series calculations near half-filling. 

The phase diagram for the t-J ladder has been dis- 
cussed by Troyer et al. Ilol. Hayward and Poilblanc 
[T^ . and Miiller and Rice |l4l|. At large J/t, the holes 
all clump together, and phase separation occurs into 
hole-rich and hole-poor regions. At intermediate J/t, 
near half-filling, a 'CISO' or Luther-Emery phase occurs, 
where the spin excitations are gapped, while there is a 
gapless charge excitation mode. Troyer et al. 0| found 
evidence of pairing between the holes in this region, to- 
gether with long-range superconducting correlations with 
modified d-wave symmetry. The spin gap is discontinu- 
ous at half-filling, as the simple magnon excitation gives 
way to a particle-hole excitation with spin. At smaller 
J/t, the phase structure appears to become more compli- 
cated, with a possible C2S2 phase appearing (two gapless 
charge modes and two gapless spin modes) T^]; while at 
extremely small J/ 1, a Nagaoka phase is expected to ap- 
pear [s^l, where each hole is surrounded by a region of 
ferromagnetic spins, forming a ferromagnetic polaron. In 
that region no spin gap occurs, and the holes repel each 
other. 

The t-Jz ladder might be expected to show similar be- 
haviour. The major difference between the models is that 
quantum spin fiuctuations are absent in the t-Jz model, 
and the system exhibits long-range antiferromagnetic or- 
der for half-filling at T = 0, whereas the t-J model does 
not. This long-range order will be destroyed at any finite 
temperature, however, and both models will then display 
similar long-range antiferromagnetic correlations. The 
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FIG. 1: The t-J ladder. 

two models should be very similar in most other aspects. 

This expectation is borne out by our numerical results. 
The phase diagram for the i-Jz ladder looks very similar 
to that of the t-J ladder, except that the critical couplings 
are about twice as large, and the Tomonaga-Luttinger 
ClSl phase extends to somewhat higher electron densi- 
ties. 

In Section ^1 we specify the model, and consider its 
behaviour in various limiting cases. In Section [llll a brief 
discussion of the DMRG method is given, and in Section 
IIVI our numerical results are presented. Our conclusions 
are given in Section Ivl 

II. THE MODEL 

The Hamiltonian of the t-Jz ladder model is (see Figure 

m 

i,a i 

-t E ^(cL.c.+i,a. + h.c.)p 

i,a,cr 

-ti_Y.P{4^„cr2.^^h.c.)P (1) 

Here the index a = 1, 2 labels the two legs of the ladder, 
i labels the rungs of the ladder, the couplings J, Jj_ are 
the strengths of the spin interactions on legs and rungs 
respectively, and t, t± are the hopping strengths on legs 
and rungs. The projection operators P forbid double oc- 
cupancy of sites as usual. A density-density interaction 
term is sometimes included as a relic of the parent Hub- 
bard model, but we do not do that here. 

In the half-filled case, with a single electron occupy- 
ing every site (n=l), the model becomes equivalent to a 
simple classical Ising antiferromagnet. The ground state 
is a doubly degenerate antiferromagnetic state (Fig. 
with energy 

i?o = -^(2J+ Ji), (2) 

where L is the number of rungs of the ladder. 

The system can be solved exactly in various limiting 
cases: 



a) 



b) 



c) 

FIG. 2: Spin configurations at 1/2 filling: a) the antiferro- 
magnetic ground state; b) an = 1 excitation; c) a domain 
wall ('soliton') excitation. 



TABLE I: Rung dimer eigenstates. 



Number 


Eigenstate 


Energy 




1 


(in>-UT))/\/2 


-J±/4 


'singlet' 


2 


TT) 1 


+Jx/4 




3 


(in) + UT))/y2 


-J±/4 


'triplet' 


4 


1 U> -1 


+Jx/4 




5 


|00) 





hole-pair singlet 


6 


(1 T 0) + 10 T))/\/2 1/2 


-ti_ 


electron-hole 


7 


(|0i) + U0))/V2 -1/2 


-t± 


bonding 


8 


(1 T 0) - |o T»/y2 1/2 


+t± 


electron-hole 


9 


(iO i) - U 0))/^/2 -1/2 


+t± 


antibonding 



1. The rung dimer limit (J/J± 0,t/J = t±/J±) 

For J_L ^ J, the system consists of independent dimers 
on the rungs of the ladder. The eigenstates on a single 
rung are listed in Table ^ The ground state is doubly 
degenerate here, with both 'singlet' and 'triplet' states 
with = having the same energy, unlike the case 
of the full t-J ladder. This degeneracy means that one 
cannot simply compute a perturbation series expansion 
about the dimer limit in this case, unless one employs de- 
generate perturbation theory, or introduces an 'artificial' 
interaction to lift the degeneracy. 
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2. The independent chain limit (J±/J —> 0,t/J — t±/J±) 

In this case we end up with two independent chains. 
The chain behaviour has been discussed by Batista and 
Ortiz (sj. They showed that the spins were Neel or- 
dered along the chain in the ground state, which could 
then be mapped onto an XXZ spin chain. Phase sepa- 
ration occurs at the (rather large) value J/t = 8 for all 
hole densities. Below that value, the system forms a gap- 
less Luttinger liquid, with correlation exponent Kp and 
charge velocity Vp. The system can be exactly solved via 
the Bethe ansatz, at quarter filling giving 



4(7r - fi) ' 



irt sin II 



(3) 



where cos^ = -J/8t. Thus for < J/8t < 1/^/2, we 
have I/V2 < Kp < 1, while for I/V2 < J/8t < 1, 
we have Kp > 1, implying that superconducting correla- 
tions dominate at large distances. At half filling, on the 
other hand, we have free fermion (metallic) behaviour 
with Kp = 1/2. 
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b) 



c) 



3. The Ising limit (t/J 0, t/J = t^/J^) 

In this limit, the model becomes equivalent to a clas- 
sical, static Ising model. Unless otherwise stated, we as- 
sume periodic boundary conditions, and an even number 
of rungs L. The ground state at half-filling is the fully 
antiferromagnetic state shown in Figure with 5*^ = 0. 
The low-energy spin excitations will consist either of lo- 
calized flipped spins ('magnons') with = ±1, as shown 
in Fig. or else domain wall ('soliton') excitations, as 
shown in Fig. analogous to the 'spinon' excitations in 
Heisenberg chains, but carrying —0 in this case. The 
total number of solitons, in the absence of holes, and as- 
suming periodic boundary conditions, has to be even for 
L even, and odd for L odd, as in the Heisenberg chain. 
The energy gaps for these two excitations are respectively 

E2b) — Eq = J + 

E2c)-Eo = J (4) 

Since there are no quantum spin fluctuation terms 
in the Hamiltonian, these spin excitations are of course 
static, in the absence of holes. Note that since the soli- 
tons carry integer spin, there is no possibility of spin- 
charge separation in this model. 

The lowest-energy charge excitation will consist of a 
single hole in the antiferromagnetic background (Fig. 
Et), with spin = ±1/2, and energy 

E3a)-Eo = ^{2J+Jj_) (5) 

The lowest eigenstates in the 2-hole sector consist of a 
bound pair on adjacent sites (Figs. |3t>ISt)i with spin 



FIG. 3: Charge excitations near to 1/2 filling; a) a 1-hole 
state; b),c) 2-hole states. 

— and energies 

E3b)-E„ = i(3J + 2Ji) 

E3c)-Eo = J(4J + Ji) (6) 

It is clear that holes will cluster together to minimize 
the number of 'broken bonds' and hence the energy, and 
the system will be phase separated in this Ising limit. 

For small but finite t/J, one can study the system via 
perturbation series calculations in t/J. Some results of 
these calculations will be shown in later sections. The 
single hole states illustrated in Fig. are localized 
states, because there is an energy barrier preventing them 
from hopping: any hop will disturb the antiferromagnetic 
alignment of the spins. At higher energy, however, there 
will be states such as that shown in Figure ^ with un- 
perturbed energy 

£^4)-£^o = J(8J+ J±), (7) 

where the hole is free to hop along the zig-zag path shown 
without incurring any further penalty in spin interaction 
energy. Assuming periodic boundary conditions, these 
states are made up of a hole plus a soliton for L even, or a 
hole alone for L odd. For finite t/J, they will form a band 
of itinerant electron states. The mobility of the electrons 
in both the chain and the ladder systems is a major point 
of distinction from the two-dimensional model, where the 
holes are 'confined', i.e. cannot move without creating a 
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FIG. 4: A hole-plus-domain-wall excitation, showing an al- 
lowed zig-zag path for hopping (solid line). 
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FIG. 5: Charge excitations for the ladder with open bound- 
aries: a) 1 hole; b) 2 holes. 



'string' of overturned spins behind them, and paying a 
penalty in spin interaction energy (s^- The states (3a) 
and (4) will mix as soon as t is turned on, and both 
configurations should be regarded as hole-soliton bound 
states. 

We will also employ open boundary conditions in our 
calculations. In this case, the holes will cluster towards 
the boundaries in the Ising limit, in order to minimize the 
number of broken bonds. In this case the ground-state 
energy is 



1 



(2./(L-l)-L./i), 



and the lowest 1-hole state (Fig. 01) has energy 

Esa) — Eq = -{J + Jj_), 

while the lowest 2-hole state (Fig. Eb) has energy 



E. 



56) 



Eo = -{2J + J^). 



(8) 



(9) 



(10) 



The single hole state is not necessarily accompanied by 
a soliton in this case. 



4- The free-fermion limit (J/t 0,t/J = t±/J±) 

In this limit, a single hole will move through the lattice 
as a free particle, and its dispersion relation is naively 
expected to be 



Ei{k) 



-2t cos k — t± 



(11) 



For simplicity, we shall restrict our remarks henceforth 
to the isotropic case t± = t, when the 1-hole energy gap 
in this limit is —3t. 

The behaviour at small J/ 1 should be similar to that of 
the t-J model. The 1-hole energy in the two-dimensional 
t-J model has been discussed by a number of authors [s^, 



|37l|38|,|39|,|40l|. At extremely smaU J/t (of order 10"^, in 
the region of the Nagaoka phase 36, 40]), the hole will 
be surrounded by a ferromagnetic 'polaron', or region of 
ferromagnetically aligned spins, with spin energy of order 
JR^ for a disc of radius R, within which the hole is fully 
mobile. In this regime, the 1-hole energy is 



E,h/t^-4 + c,{jy/' 



(12) 



where ci is a numerical constant. At larger J/t, how- 
ever, the lowest-energ^y states correspond to Brinkman- 
Rice 'string' states 37, "38], where the hole is confined by 
a 'string' of overturned spins within an antiferromagnetic 
background. This picture gives a 1-hole energy 



Eu^/t^-a, + b^{Jf/^ 



(13) 



where ai < 4. Shraiman and Siggia '38l| estimated ai = 
2\/3 and bi — 2.74 for the t-Jz model, quite close to the 
values found numerically by Barnes et al js^ , and White 
and Affleck "i^. 

For a one-dimensional chain, no Nagaoka phase occurs, 
but for a ladder the condition is once again met that 
holes can ho p ar ound closed loops, and a Nagaoka phase 
is expected [l3. White and Affleck |40| have pointed 
out that the ferromagnetically aligned spins inside the 
polaron will orient themselves in the xy plane, rather 
than along the z axis, because this costs an exchange 
energy of J/4 per bond, rather than J/2. For our ladder 
model, the ferromagnetic 'polaron' region will be one- 
dimensional. If the polaron region spans L rungs of the 
ladder, the cost in magnetic energy will be 



Em — -JL 



(14) 



Assuming the 1-hole wavefunction to vanish at the edges 
of the polaron, it will take the form 



Hk) =Aj2cos{kn)i\n,l) + \n,2)), k = 



L 



(15) 



where \n,a) denotes a hole at rung n on leg a of the 
ladder, and A is a normalization constant. The kinetic 
or hopping energy is then easily found to be 



E 



K 



-'it + t{-f 



Minimizing the total energy E^^ 
to L we find 



{L^od) (16) 
Em+Ex with respect 



Eih/t 



-3+J(^)^/^ 
-3 + 3.35(J/i)2/3 



(17) 
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For intermediate J/t, the Brinkman-Rice string picture 
|37l l38j| will give a qualitatively similar behaviour 



(18) 



where 02 < 3. It is doubtful whether approximations 
such as those used by Shraiman and Siggia .^^l to cal- 
culate the coefficients 02 and 62 for the 2D model are at 
all accurate for the ladder, and we will not attempt an 
explicit calculation of these coefficients. 



using (|21|l . where n is the electron density. This formula 
applies for either periodic or open boundary conditions, 
in principle. The velocity may be measured indepen- 
dently using the excitation energy of the lowest state of 
momentum 27r/L for periodic BC 



Einm = l)-Eo = 



(23) 



For open BC, the momentum spacing is n/L, and the 
corresponding formula is 



A. The Effective Hamiltonian 

In the regime of physical interest, the t-J ladder is 
believed to be in a CISO or Luther-Emery phase, with 
gapped spin excitations and gapless charge excitations 
cor resp onding to bound hole pairs. Several authors [inL 
[21I [2^ have discussed an effective Hamiltonian to 
describe these bosonic excitations, which would capture 
the low-energy physics of the model. Here we merely 
summarize their results. 

In the recent analysis of White, Affieck and Scalapino 
[2^ . they use a bosonization technique to construct the 
low-energy effective Hamiltonian 



H - fiL 



dx[KpIll 



1 



Kn dx 



(19) 



where is a chemical potential, Hp is the momentum den- 
sity conjugate to 9p, Vp is the velocity of the correspond- 
ing gapless low-energy excitations and the parameter Kp 
controls the correlation exponents. The two parameters 
Vp and Kp must be extracted from numerical data. 

White et al. '2^ show that the finite-size scaling be- 
haviour of a general low-energy excitation is 



(24) 



Hence one can obtain estimates for Kp and Vp separately. 

One may also use 'twisted' boundary conditions with 
the wave function acquiring a phase $ at the boundary, 
corresponding to an Aharonov-Bohm flux threading the 
one-dimensional ring formed by joining the ends of the 
ladder together. This increases the ground-state energy 
by 



Eq Eq 



SnVpKp^^ 



(25) 



and allows one to determine the combination VpKp. This 
approach was used by Hayward and Poilblanc (12] . 

Finally, White et al. [23| discuss how one may es- 
timate Kp from the decay of Friedel oscillations in the 
system. Friedel oscillations are density oscillations pro- 
duced near the boundary of an open ladder, which decay 
with a power law corresponding to the exponent Kp with 
distance into the ladder. Using their bosonization analy- 
sis and a conformal transformation. White et al [2^ show 
that for a system of length L the density at site (rung) j 
should vary as 



E-Eq = -2pn- 



2ttv 



'' ^Kpfn^ + -J— + y] k{nLk + nRk)] 

4: jV o 



(20) 

where Eq is the ground-state energy for a given density 
n, riLk and nuk are occupation numbers for left and right 
moving states of momentum 27rfc/L, L is the number of 
rungs of the ladder, and p and m are integer- valued quan- 
tum numbers. The total charge relative to the ground 
state is Q = —2p, and the other quantum number m 
measures the "chiral charge" . 

By measuring the ground-state energies for three dif- 
ferent charge states with AQ = ±2, one can determine 
the ratio Vp/ Kp: 



E{p ^ 1) + E{p ^ -1) - 2Ea 



KpL 



(21) 



This is directly related to the electron compressibility k 
of the two-leg ladder, generally defined as 



1 



1 d^E 



2L dn^ 2K„ 



(22) 



Ccos{2imj + (3) 
[(2L/7r) sin(^j7L)]^p 



(26) 



where C, /? are constants, and 27rn is the minimum wave 
vector of Friedel oscillations in the CISO phase, corre- 
sponding to two holes per wavelength. This is the same 
wave vector that would occur for a one-component spin- 
less hard core bose gas, made up of tightly bound hole 
pairs. 

Schulz has argued in the case of the t-J ladder that 
Kp should take the universal value Kp = 1 at half-filling; 
but his argument starting from the rung dimer limit is 
not necessarily applicable to the t-J^ model. 

Note that the picture may be complicated by the oc- 
currence of a gapped charge density wave (CDW) phase 
at special commensurate filling factors such as n = 0.5 
or n = 0.75. In this case, the charge density oscillations 
persist at all distances from the walls, corresponding to 
broken translational symmetry and a regular pattern of 
hole placements. This possibilit y h as been discussed by 
Riera et al [33 and White et al 
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block size by a whole rung at each DMRG step as in Fig. 
El We then perform a finite lattice sweep at the fixed lat- 
tice size, using the previously developed system blocks to 
improve the accuracy of the calculation. Once a target 
size has been completed with a finite lattice sweep, the 
infinite lattice algorithm is resumed to proceed to a larger 
target lattice size whereupon the sweeping procedure is 
again employed. This way we obtain accurate finite lat- 
tice method results for a number of lattice sizes, e.g., 
iV = 4, 8, 12, 16, 24, 32, 64, 128, 256 lattice sites. Wher- 
ever necessary, the data are extrapolated to the bulk limit 
A'^ — s- oo by fitting a low-order polynomial in to the 
finite-lattice data. 



IV. RESULTS 



O 



o 



o 



o 



FIG. 6: Infinite lattice construction of the two- leg ladder. 
Two whole rungs (four sites — open circles) are added per 
DMRG step. The system and environment blocks increase 
by one rung (two sites) per DMRG step. 



III. METHOD 

The model is solved using the Density Matrix 
Renormalisation Group (DMRG) method 42] . Cal- 
culations have been performed at both half-filling and 
quarter-filling. Calculations have been performed with 
both periodic and open boundary conditions in the hor- 
izontal direction, but most of the results use OBCs be- 
cause the convergence is much better in this case. 

In the half-filled case we have calculated the ground 
state as well as one, two and four-hole excited states. In 
the one-hole case we have calculated S*^ = 1/2, 3/2 and 
5/2 states. We have calculated energies of these states as 
well as density profiles {fii). The infinite lattice algorithm 
|4^ is employed in this case, using typically M = 60 
states per block and symmetry sector, which corresponds 
to a total of about 300 system block states. The lattice 
build-up phase is illustrated in Fig. |^ By running tests 
for a number of different values of M, we established 
that the energies of these states have been resolved with 
sufficiently high accuracy for our purposes. 

In the quarter-filled case, it is necessary to use the fi- 
nite lattice method [i^l in order to achieve reasonably 
converged energies. We use the infinite lattice method 
to build a ladder to a given size, increasing the system 



We have calculated numerical DMRG results for this 
model for the isotropic case, J = J± = Jz,t = t±, and 
for various Jz/t, on lattices of L rungs with L even. Our 
discussion henceforth will be limited to this case. 



A. Near Half-Filling 

We begin by considering states with a finite number of 
holes doped into the half-filled system (n=l). 



1. Single-hole states 

Figure [7| shows DMRG estimates of the energy gap for 
a single hole with spin — 1/2, as a function of Jz/t, for 
both periodic and open boundary conditions. The results 
of a series calculation are also shown for the periodic case. 
The series estimates were obtained assuming the Nagaoka 
form 



fit) = Eih + 3i - 5^1/3 



t^O. 



(27) 



Accordingly, the series in t for f(t) was Euler transformed 
z — t/(l + t) to bring the singularity to z = 1; then a 
further change of variable was made 



1 



(1 



\l/3 



(28) 



so that one may treat the function {l — 6)f{5) as analytic 
in 6 near d — 1; and finally, differential approximants 
were used to extrapolate this function to 5 = 1. It can 
be seen that the series results agree very well with the 
DMRG results for J — z/t > 0.3, but run a little lower 
below that. 

It can also be seen that there is a difference between 
the DMRG results with periodic and open boundary con- 
ditions. As discussed in Section III 31 this is because 
the state with periodic boundary conditions is actually a 
hole-soliton bound state. The difference between the two 
energies is shown in the inset to Figure 
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— Series expansion 
• -• DMRG PBC M=80 N=64 

DMRG OBC M=80 N=64 

DMRG OBC with end term, IVI=60, N=132 




0.15 




0.05 



FIG. 7; (Colour online) The energy gap E^h/t = (-Bi - Eo)/t 
for one hole, as a function of Jz/t- Filled circles, open bound- 
ary conditions; open squares, periodic boundary conditions; 
dashed line, estimates from Ising series expansion. 




FIG. 8: (Colour online) Low energy spectrum of the t-J (left) 
and t-Jz (right) model on a 4-rung ladder with 7 electrons. 



At smaller Jz/t, one is tempted to look for a transi- 
tion to a Nagaoka phase [s^, i.e. a phase with ferro- 
magnetic spin order or at least a ferromagnetic polaron 
bubble around the hole. In the t-Jz model, however, the 
situation is rather different from the t-J model, where the 
full spin-rotation symmetry applies. The scenario is il- 
lustrated in Figure|Sl where we show the complete 1-hole 
energy spectrum of the t-J and t-Jz models on a 4-rung 
cluster. The left-hand panel shows the t-J case, where 
at J = the ground state has maximal spin S = 7/2 
(ferromagnetic). As J increases, at some critical cou- 
pling Jc this state crosses another state with minimal 
spin S — 1/2, which becomes the ground state there- 
after. The right-hand panel shows the t-Jz case. At 
Jz — 0, the system is rotation symmetric, the Nagaoka 
theorem applies and the ground state again has the max- 
imal possible spin S'max and degeneracy 2S'inax + 1- As 
soon as Jz becomes finite, however, the symmetry is bro- 
ken, the multiplet splits into its components, and the 
state with minimal S*^ = ±1/2 becomes the ground state. 
Clearly, there is no level crossover in the t-Jz case. All 



FIG. 9: The 1-hole bandwidth W/t as a function of Jz/t, 
estimated from the Ising series expansion. 



we can expect for increasing Jz is a continuous fading-out 
of the ferromagnetic correlations. This accords with the 
smooth approach of the 1-hole energy to — 3 as Jj ^ 
in Figure 

A comparison of the numerical data at small Jz/t with 
the theoretical prediction H17|) gives the following results. 
A Dlog Fade analysis of the series gives a rather inaccu- 
rate estimate of the exponent in the range 0.5 - 0.7, which 
would be consistent with 2/3. A fit to the DMRG data 
over the range — 0.1, however, gives 



Eih 



= -3 + 1.444 



0.53 



(29) 



with an exponent much closer to 1/2 than 2/3. This 
is illustrated by the difference between the DMRG data 
and the series extrapolation in Figure [7| at small Jz/t. 
Thus the data do not seem to behave in accordance with 
theoretical expectations in this instance. 

Figureinishows the bandwidth for the 1-hole state with 
= 1/2 predicted by the Ising series expansion, as a 
function of Jz/t. We sec that the predicted bandwidth 
rises as Jz/t decreases, as one might expect, until at 
Jz/t ~ 0.3 it reaches a peak, and begins to drop to- 
wards zero. Similar behaviour in the region of small J/t 
has been predicted for the 2D t-J model by Martinez 
and Horsch and Liu and Manousakis and con- 
firmed by series calculations The theoretical pre- 
dictions assume long-range antiferromagnetic order, and 
construct an effective Hamiltonian in which the spin dy- 
namics are treated in linear spin-wave theory, and the 
holes are treated as spinless fermions, following Schmitt- 
Rink, Varma and Ruckenstein [i^. The results are con- 
sistent with the string picture of Brinkman and Rice '3 7*1 , 
and appear to show a sequence of narrow 'string' exci- 
tation peaks in the spectral function, corresponding to 
holes bound in a linear potential. The effects are even 
more marked in the t-Jz model. 

It would be very interesting to check the actual be- 
haviour of the bandwidth by means of DMRG estimates. 



FIG. 10: (Colour online) Left: Two-hole binding energies as 
a function of Jz/i, calculated with OBC and finite boundary 
potentials. Filled circles: = 0: filled squares: 5"^ = 1. 
Right: Average electron number per rung as a function of 
distance along the ladder, for 2-hole states on a lattice of 34 
rungs: a) J^/t — 0.5; b) J^/t = 1.0. The lines are fits using 
equation 



012345678 
J,/t 

FIG. 11: (Colour online) Four-hole binding energy as a func- 
tion of Jz/t, calculated with OBC and finite boundary poten- 
tials. Insets: Electron density profiles at characteristic values 
of Jz/t, together with fits to the corresponding Friedel expres- 
sion. 



but unfortunately our codes cannot distinguish different 
momentum eigenstates. 



2. Two-hole states 

In studying 2- and 4-hole states and looking for the 
phase separation boundary, we wanted to avoid the holes 
sticking to the boundary and thus obscuring the bulk 
physics. Accordingly, we have added an extra 'boundary 
potential' term at the sites adjacent to the boundary: 



Hf 



4 ^4,a(^hT- 

i,(T 



(30) 



As it turns out, this makes a negligible difference to 
the energy gaps in the region of interest. 

In the 2-hole sector, the first quantity of interest is the 
2-hole binding energy 



Ef 



2Ei — E2 — Eq 



(31) 



Figure [TUl shows the binding energy as a function of Jz/t 
for both 5^ = and — 1, computed with open bound- 
ary conditions. It can be seen that a 2-hole bound state 
with = forms for large Jz/t, but the binding energy 
vanishes at Jz/t ~ 0.60, and below this point the holes do 
not bind. In the — 1 channel, there is a small binding 
energy at large Jz/t, which vanishes at Jz/t ~ 1.8. 

Similar results have been found for other members of 
this fami ly o f models. For the 2D t-Jz model, Riera and 
Dagotto [33 estimated that the holes become unbound 
at Jz/t ~ 0.18; while for the t-J ladder, Jurecka and 
Brenig [s^l suggest that hole binding vanishes at J/t ~ 
0.5. There is binding in the S' = 1 channel for the t-J 



ladder, but it is considerably smaller than in the 5 = 
channel. 

Wherever the 2-hole binding energy is positive, we ex- 
pect a continuous band of 2-hole bound states to appear 
above the ground state; so that within the 2-hole sector, 
particle-hole excitations ('excitons') can occur down to 
zero energy. In other words, excitons with ~ are 
gapless, as found by Troyer et al 10] for the t-J ladder. 

We note that the spin energy gap between the — 1 
and — 2-hole states also drops to zero at Jz/t = 0.6. 
It always remains below the — 1 'magnon' energy 
3 J/2; or in other words, an exciton with S*^ = 1 in the 2- 
hole sector has an energy substantially below that of the 
magnon. Therefore the spin gap drops discontinuously 
as one moves away from half-filling, which again agrees 
with the behaviour found by Troyer et al. 10] for the t-J 
ladder. 

The right panels of Figure [TUl show profiles of the elec- 
tron density in the 2-hole state as a function of distance 
along the ladder, calculated for a lattice of 34 rungs with 
open boundary conditions. At low Jz/t, the two holes 
are unbound and separated from each other, forming 
two distinct troughs in the density profile (Fig. ITUI) . At 
Jz/t ~ 0.6 a sharp crossover takes place, where the two 
holes bind together and form a single trough in the den- 
sity profile. As it happens, both profiles can be fitted 
quite well with a Friedel oscillation form 



Co 



ci sin[27r(n+ 1/2)j7(£+ 1)] 
[{L + 1) sin[^j7(L + l)]]^p 



(32) 



where n = 2 and Kp = 0.87 in case a) and n = 1 and 
Kp = 0.92 in case b). 
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FIG. 12: (Colour online) Four-hole binding energy as a func- 
tion of XjN = 1/{2L), at J^/t = 1, 2 and 3. 

3. Four-hole states 

The 4-hole binding energy for 5^ = is shown as a 
function of Jz/t in Figure ITTl It is finite at large Jz/t, 
and vanishes smoothly at Jj/f ~ 4. This is taken to mark 
the point where phase separation occurs. The smooth be- 
haviour of the binding energy may indicate a Kosterlitz- 
Thouless transition at this point. The insets show exam- 
ples of the electron density profiles as functions of dis- 
tance along the ladder for the 4-hole ground state. For 
Jz/t = 0.5, we see four separated holes; for Jz/t = 1.0, 
we see two separated hole-pairs; and for Jz/t = 6.0, we 
see a single 4-hole cluster, reinforcing the picture given 
above. Note that the Friedel form (|32(l again fits the 
density profiles remarkably well, upon choosing the ap- 
propriate value for n. The fits in the pair binding region 
indicate that the parameter Kp lies approximately in the 
range 0.8 - 1.0 near half-filling. 

Next, we attempt to measure the isothermal compress- 
ibility in the non-phase-separated region, whose inverse 
is 

= n^N[Ei + Eo- 2E2\/A (33) 

where = 2L is the number of sites, and n — {N ~2)/N 
the electron density. Figure IT^ shows K~^n~^ as a func- 
tion of 1/A^, for some couplings Jz/t < 4, fitted by a 
quadratic polynomial in 1/iV. The extrapolations to the 
bulk limit N oo are compatible with zero, indicating 
that the inverse compressibility vanishes at half-filling. 
This seems surprising at first sight, but the same re- 
sult applies for the ID t-J chain at the integrable point 
[47|. Siller et al. also found a tiny result for the 

inverse compressibility of the isotropic t-J ladder near 
half-filling. The result can be understood as follows. The 
dispersion relation in the conduction band flattens out 
near the band edge, and hence the dispersion relation 
of a single hole at half filling is quadratic in momentum 
g, as illustrated by the free fermion equation 1111 or the 
calculation of Sushkov for the t-J ladder. The dis- 
persion relation for a bound hole-pair bosonic excitation 
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FIG. 13: The odd-even gap [E{nh = L + 1,5*^ = -1/2) + 
E{nh = L - 1,5" = -1/2) - 2E{nh ^ L,S^ = 0)]/2 at 
quarter-filling, as a function of Jz/t. 

will likewise be quadratic in q. The vanishing of any term 
linear in q implies that Vp ^ and Eq. (|22|l then leads 
to K-^n-'^ = 0. 



B. Quarter Filling 

Next we consider states at or near the quarter-filling 
point {n = 1/2). Reliable results are very hard to obtain 
for Jz/t > 3. Phase separation presumably occurs in this 
region, so that the ground state becomes inhomogeneous, 
and very difficult to treat by DMRG methods. For the 
most part, our results will be for smaller Jz/t. 

1. Odd-even gap 

Figure 1131 shows the energy gap between states with 
odd and even numbers of holes, [E{nh = L -\- \,S^' — 
-1/2) + E{nh = L-1,S^ = -1/2) - 2Einh = L,S' ^ 
0)]/2, where Uh is the number of holes, as a function 
of Jz/t. It can be seen that for Jz/t less than about 
2.0, this gap is small and positive, but is consistent with 
zero beyond that point, within substantial errors. This 
behaviour is very similar to the two-hole gap shown in 
Figure El In fact it appears that the energy surface is 
quite smooth in this region, with little if any 'corrugation' 
between states with even or odd numbers of holes. This 
indicates that pairing effects between holes are weak or 
nonexistent in this region. 

2. Spin Gap 

Figure [TH shows the spin gap [E{nii ^ L,S^ = 1) — 
E{nh = L,S^ = 0)] as a function of Jz/t. The insert 
shows examples of the finite-size scaling behaviour of this 
quantity as a function of 1/A^, and the linear extrapola- 
tions made to the bulk limit. It can be seen that the 
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FIG. 14: The spin gap [E{nh = L, 5^ = 1) - E{nh = 1,3" ^ 
0)] at quarter-filling, as a function of Jz/t. 




1/N 



FIG. 15: The two-hole gap [E{nh = L -f- 2, 5"^ = 0) + E{nh = 
L - 2, 5^ = 0) - 2E{nh = L,S'^ ^ 0)]/2 at quarter-filling, as 
a function of 1/N = 1/(2L), at selected values of Jz/t. 



results are consistent with zero over the whole region 
Jz/t < 4. Convergence for J^/t > 2 is not very good. 
Thus for Jz/t > 2.3, we appear to be in the ClSl phase 
discussed by Hayward and Poilblanc 12] for the t-J lad- 
der, where both the charge gap and the spin gap vanish. 
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FIG. 16: The estimated two-hole gap in the bulk limit 
[E{nh + L + 2,S'' =0)+ E{nh = L - 2, 5" = 0) - 2E{nh = 
L, = 0)]/2 at quarter-filling, as a function of Jz/t. 

in the electron density profiles and correlation functions. 
The density profile on an odd-rung lattice at quarter fill- 
ing shows some evidence of rung-to-rung oscillations, but 
they die away as iV — > oo, and so do not seem to corre- 
spond to the CDW. 

Examples of the density-density correlation functions 
are shown in Figure [T7I with one site fixed at the middle 
of the ladder, and the distance to the other site varying, 
including both sites on the same leg and those on the 
opposite leg of the ladder. There is a positive correlation 
for nearby sites, and a Friedel rung-to-rung oscillation 
near the boundary. In the intermediate region, there is 
no strong evidence of either a rung-to-rung alternation or 
a checkerboard alternation. One major qualitative differ- 
ence is obvious, however, between Figures ITTk) andlTTh). 
At Jz/t = 2, the correlation function is virtually zero at 
intermediate and large distances, as one would expect of 
a solid phase: this is consistent with the CDW behaviour. 
For Jz/t = 3, on the other hand, the correlation function 
has a substantial dip at intermediate distances, which 
actually increases with increasing lattice size. This be- 
haviour is more reminiscent of a liquid or a gas. 



3. Two-hole gap 

Figure [TKl shows the two-hole energy gap [E{nfi = L + 
2, S*^ = 0) + E{nh = L-2, = 0)-2£;(?i,, = L,S^ = 0)] 
as a function of 1/iV at selected values of Jz/t. For Jz/t < 
2.3, the gap appears to scale to a finite value in the bulk 
limit N oo, whereas the limit is consistent with zero 
from there on. Figure ITHl shows the estimated bulk limit 
as a function of Jz/t. The finite charge gap for Jz/t < 2.3 
corresponds to a cusp in the energy density versus filling 
factor surface, and is evidence of the formation of a CDW 
state at this commensurate filling factor, according to 
White et al. 

We have looked for further evidence of the CDW state 



V. DISCUSSION 

As expected, the phase behaviour of the t-Jz ladder 
at T = appears to be broadly similar to its counter- 
part, the t-J ladder. A schematic diagram is shown in 
Figure^! At half- filling, the system is equivalent to an 
Ising antiferromagnet, with long range magnetic order; 
but at any finite temperature or hole density, the order 
parameter will vanish, leaving only finite-range antifer- 
romagnetic correlations. 

Near half-filling, phase separation into hole-rich and 
hole-poor regions occurs beyond Jz/t ~ 4.0, which is 
about half the critical value for the t-Jz chainJSJ, but 
twice the critical value for the full t-J ladder [l3|. Below 
that, there is a CISO or Luther-Emery phase, where the 
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FIG. 17: The density-density correlation function as a func- 
tion of distance, with one site fixed at the middle of the ladder 
and the other varying, at quarter filling. Triangles: L = 22; 
squares: L — 88. Case a): Jz/t = 2.0; b): J^/t = 2.6; c): 
Jz/t = 3.0 



holes bind into pairs and the charge gap vanishes, but the 
spin gap remains finite. The spin gap is discontinuous 
at half-fiUing, because the exciton spin gap is less than 
the magnon gap. These features are identical with those 
found for the t-J model We have not explored the 
question whether there are superconducting correlations 
in this region, but they are expected to be present if 
Kp > 1/2. 

The binding between hole pairs vanishes in the — 1 
channel below J^/t ~ 1.8, and in the = channel 
below Jz/t ~ 0.6. Below that, the holes repel each other, 
and we presumably enter the C2S2 gapless phase dis- 
cussed by Miiller and Rice 0| for the t-J case (although 
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FIG. 18: Schematic phase diagram for the t-Jz ladder. 



exact diagonalization data up tp iV = 14 are somewhat 
equivocal as to whether the = 1 states are degener- 
ate). These different regimes are clearly evident in the 
hole density profiles on a finite lattice. 

In the Nagaoka limit as Jz/t ^ 0, a finite-lattice calcu- 
lation reveals rather different behaviour in the two mod- 
els in the one-hole sector. In the t-J case, a state with 
maximal spin S crosses over the 5 = 1/2 state to be- 
come the ground state as J/t — > 0. For the t-J^ sys- 
tem, on the other hand, the — 1/2 state remains 
the ground state at all Jz/t, and the Nagaoka limit is 
approached smoothly. States with higher become de- 
generate with the — 1/2 state in the limit, to form a 
multiplet with maximum total spin S. Any ferromagnetic 
correlations in the ground state appear to be transverse 
to the z direction, as arg ued by White and Affleck 4(|. 
Our DMRG results for the one-hole gap appear to de- 
crease as (Jj/t)^/^, rather than {Jz/tY^^ as predicted by 
theory; the theory needs re-examination for this case. 

Some crude estimates of the effective Hamiltonian pa- 
rameters in the Luther-Emery phase were made at half- 
filling. Friedel-type fits to the density profiles give Kp ~ 
0.8 — 1.0, although it is not clear that the Schulz argu- 
ment that Kp = 1 at half-filling is applicable to the 
i-Jz case. We found that Up = in that limit, a general 
result which should apply to chains and ladders for both 
models. 

We have also made some studies of the quarter-filled 
case. For Jz/t < 2.3, the system appears to be in a CDW 
state, with a finite two-hole gap and solid-like density- 
density correlation functions, although we could not dis- 
cover any particular pattern of charge density waves. 
This behaviour was predicted for a commensurate fill- 
ing factor by White, Affleck and Scalapino [23|. We find 
that the spin gap vanishes in this phase, and in fact is 
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compatible with zero at all couplings. For J^/t > 2.3, 
any pairing effects are weak or nonexistent, and the sys- 
tem appears to be in a ClSl Tomonaga-Luttinger phase 
with a single gapless spin and single gapless charge mode, 
similar to that discussed for the t-J ladder at densities at 
or below quarter-filling by Hayward and Poilblanc [T^ . 
The calculations become very unstable beyond J^/t 3, 
and we take this as indicating the onset of phase separa- 
tion again at larger Jz/t. At very small electron densities, 
an electron-paired phase should exist. These findings are 
indicated schematically in Figure ITSl 
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